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We investigate whether or not Kondo insulators undergo 
a magnetic field induced metal-insulator transition in one di- 
mension at half filling using both a density matrix formulation 
of the numerical renormahzation group and a new analytical 
approach. Contrary to expectations, the quasiparticle gap 
never vanishes at any field and no metal-insulator transition 
is found. We discuss generalizing our result to three dimen- 
sions and to the asymmetric Anderson lattice. Our result is 
consistent with recent experiments. 



In a strong magnetic field an ordinary narrow gap 
semiconductor becomes metallic when the field is com- 
parable to the size of the gap. Is this also true for the 
class of rare earth compounds known as Kondo insula- 
tors, whose insulating gap is due to interactions between 
conduction electrons and localized f-electrons Sim- 
ple arguments j^-^ suggest that such a metal-insulator 
transition (MIT) should occur when the applied mag- 
netic field h exceeds the exchange coupling J but is less 
than the conduction electron bandwidth. In this case, 
one expects the f-spins to be completely polarized and 
decoupled from the conduction electrons, leaving an in- 
completely polarized band of metallic conduction elec- 
trons. This scenario is supported by mean field calcula- 
tions such as those with slave bosons However, here 
we show that such expectations are not correct when spin 
fluctuations are taken into account. 

We have studied the behavior of the insulating gap in 
a one-dimensional Kondo insulator as a function of mag- 
netic field. Using the density matrix formulation of the 
numerical renormahzation group @,| as well as a new 
analytic approach, we find that the gap, while greatly 
reduced by the applied magnetic field, never vanishes at 
any value of h. This unexpected result is due to spin 
flip scattering across the Fermi surface which removes 
degeneracies and keeps the gap open. We argue that this 
result also holds in three dimensions and for the asym- 
metric Anderson lattice where there is no particle-hole 
symmetry. Our results are consistent with recent resis- 
tivity and Hall measurements on Ce3Bi4Pt3 in pulsed 
flelds up to 61 T which do not see the signatures ex- 
pected for a transition from a semiconductor to a good 
metal This material acts like a semimetal or a dirty 
metal at the lowest temperatures (T < 5 K) and high- 
est fields {h > 50 T), even though the field exceeds the 
experimentally deduced value of the gap (A ~ 37 K) 



The ID Kondo lattice has spin- 1/2 conduction elec- 
trons that hop from site to site with an on-site spin 
exchange between a local f-electron and the conduction 
electron on that site. Thus the Hamiltonian is 



iJ = i^cl^a+ia + h.c. 

i.<7 

+ J2 JS.f ■s.c-hY^ {SL + S!f) , (1) 

i i 

where the conduction electron spin density on site i is 



Sic = 4a(^/2)a/3Ci/3, the f-electron spin density is 

Sif = J2af3 fLi^/'^)af3fil3^ ^afi are Pauh matrices and h 
is the external magnetic field, chosen along the z-axis. 
We set the hopping matrix element t = 1 and choose 
J to favor antiferromagnetic ordering (J >0). To study 
the Kondo insulator, we restrict ourselves to half filling 
where the total number of conduction electrons N equals 
the number of sites L. 

We consider a system to be insulating when the quasi- 
particle gap Aqp = /iAT+i — 7^ 0. Here ^ln = 
En — En-1 is the chemical potential and En is the 
ground state energy with N electrons. We have calculated 
Aqp using both the density matrix renormahzation group 
(DMRG) algorithm P,^|Jl0|] and a new analytic technique. 
We first describe the results of our numerical work. 

The DMRG approach is a real space technique which 
has proven to be remarkably accurate for the Kondo lat- 
tice . We used the finite system method with open 
boundary conditions in which there is no hopping past 
the ends of the chain. We studied lattices of size L—6, 8, 
16, and 32, keeping up to 100 states. The energies were 
extremely accurate for J ^ t, with typical truncation er- 
rors of order 10~^ for J=10. For J^t, the f-spin degrees 
of freedom lead to a large number of nearly degenerate 
energy levels. As a result, the accuracy was significantly 
reduced, with truncation errors of order 10^^ for J=0.75. 

To determine A^p, we calculated the ground state en- 
ergies with N = L and iV = L -I- 1 conduction electrons 
pl| for all possible values of the total with ft, = 0. 
Since the magnetic field term in the Hamiltonian com- 
mutes with the Hamiltonian, we can find the energies in 
the presence of a field by adding —hS^^^ to the zero field 
energies. We then use the lowest energies (i?]^™(S'fot) 
and i?]^™ (5fot)) as a function of Sf^t to calculate A^p. 
In the presence of a field, finite size effects do not limit 
the size of the gap because the field can tune the levels 
En and En+i to be infinitesimally close. For example, 
for free electrons on an 8 site lattice, the minimum Aoo 
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is zero to within the roundoff error of 10~^^. In Fig. 1 
we show the quasiparticle gap as a function of h. The 
sawtooth osciUations occur as the ground states, En and 
En+1, change their polarizations with increasing h. 

Notice that at high fields, the gap is proportional to 
h. This occurs when the spins of the half- filled lat- 
tice are completely polarized by the field. To calcu- 
late i?Ar+i, we add an electron whose spin is opposite 
to the field. The energy required to do this increases 
linearly in h, and hence Aqp is proportional to h. This 
can also be understood in terms of splitting the up(+) 
and down(— ) spin conduction electron bands. For fully 
polarized f-spins, the conduction electron energies are 
E± = -2tcos{k) ± (J/4 - h/2). The two bands are 
completely separated at a field hs = J/2 + At and so, 
for h > hg, Aqp ^ h. Our expression for hg is only an 
estimate since, as we shall see, the f-spins are not fully 
polarized until the bands separate. 

For a given value of J, we can find the minimum quasi- 
particle gap. A™™, as a function of h from plots similar 
to Fig. 1. Roughly speaking, the minimum occurs at the 
largest field in which the up and down spin bands are not 
completely separate. This field, hmm, is of order the con- 
duction electron bandwidth, i.e., hmin ^ hs ^ 4i. Notice 
that the minimum does not occur at h J or at h Aqp 
as one might expect. Fig. ^ shows A™" as a function of 
J. Fitting the curves for J < 2 to A™" ~ J^ we find 
that b increases from 2.05 for L=6 to 2.16 for L=32. The 
power law form indicates that the gap does not vanish at 
any finite value of J. 

To see why the gap remains open, note that the f-spins 
are not completely polarized for h < hg. This allows spin 
flip scattering across the Fermi surface between the con- 
duction electrons and the f-spins. It removes degenera- 
cies at the Fermi energy and prevents the quasiparticle 
gap from closing. The inset of Fig. |^ shows the con- 
duction band structure for a nonzero h. Note that the 



Fermi points satisfy k 



R+ 



TT and kp 



kp^ 



for h < hg. Thus scattering processes with momentum 
transfer 2kp = tt will be important, {kp is the Fermi 
wavevector in zero field.) 

We have verified this scenario using a new analytic ap- 
proach in which we keep two types of scattering processes: 
(1) those with no momentum transfer and no spin fiip 
and (2) those with spin flip scattering with tt momentum 
transfer. We call these SI and S2, respectively. These 
processes are not restricted to the Fermi surface. 

We use the Popov-Fedotov functional integral ap- 
proach to enforce the constraint of 1 f-electron per site 
Ip^ . Our strategy is to start from the partition function 
exp(— /3iJ) and integrate out the f degrees of freedom 
self-consistently to obtain an effective action for the con- 
duction electrons. This yields the conduction electron 
propagator whose poles determine the T=0 excitation 
spectrum. The gap is the smallest excitation energy. We 



find a nonzero gap for all values of the field. 

We adopt a basis of Right and Left variables for the 
conduction and f-states. If ■0 is a c or / operator, then 



IpRaik^Lo) 

V'La(fc,a;) 



ipa{k, uj) k > 
otherwise 

1pa{k ~ TT^Uj) k > 

otherwise 



(2a) 
(2b) 



Thus, 2kp momentum scattering processes link R and 
L variables, while zero momentum terms involve R or L 
states only. We write the partition function, in a Fourier 
representation, as follows: 



Z 

Soc 
Sof 
Sfc 



VcDcDfVfe' 



UJ 

-y 



dk 

2^ 



c{k, uj)Doc{k, oj) 



rlk - 

— f{k,cu)Fof{k,u;) 



(3a) 
(3b) 

(3c) 



1 



dkidk2 



/(fcl,C.l)Cijf/(fc2,C.2) (3d) 



/32 ^ 7o (27r) 



where c and / are now vectors of Grassman variables of 
the form tp = {ipR+^ V'l+i V'fi-j '4'l-)- w is a fermion Mat- 
subara frequency, i.e., ui — 27r(n -I- 1/2)//?. Afc = ki~k2, 
Auj = uji — UJ2, and (3 is the inverse temperature. Soc 
and Sof are the actions for the free conduction electrons 
and the f-states respectively, and D^^ and F^^ are the 
corresponding (imaginary time) propagators: 



Do = 
Fo = 



ioj — h — 2t cos(A:)(T^ 

iu! + h — 2tcos{k)a^ 

iLU - m/2P -h 

iuj - iTr/2/3 + h 



TT Sfc contains the coupling between the conduction elec- 



trons and the f-states: 



Ak 



A = 



A 2B__ 
2B. -A 



Baf) = 



i?4 



Rap.Q RaP,TT 
RajJ.TT RaPfi 

R ,0 R++..Tr ^ R ,7r 

R ,7r R++fi ~ R ,0 



(4) 



Doj Fo and are 4x4 matrices while A and Bap are 
2x2 matrices. Rap.o represents scattering processes with- 
out spin flip ("0"): 

7-> J dk^dkji , , . . 

= 7-5 > / 1^ .2 {CRa{k3,L^)CB,0[k4,UJ + Auj) 
4/3 Jo {2ny 

+ CLaik3,Uj)cLl3ik4,UJ + ALu)}S{ki - fcs - Afc) (5) 

and Rap.-n describe spin flip scattering across the Fermi 
surface (tt): 
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^a/3,7r = |^2^ J -j^;^{CRa[k3,UJ)CL0[ki,UJ + Auj) 

+CLa(fc3, t^)cRi3{ki, w + Aa;)}^(A:4 - ks ~ Ak - tt) (6) 

where S enforces momentum conservation mod27r, i.e., 
umklapp scattering is included. 

Since we are interested in the c fields, we integrate out 
the / fields. This leaves an effective action Seff for the 
conduction electrons: 



Z 



g «e// ^ I g-s„c+in(det(j=^„+c)) 



An expansion of the ln(det..)) around Fo would corre- 
spond to perturbation theory in terms of the bare 
propagator. This is not a good starting point since it 
ignores the effect of the conduction electrons on the f- 
states. Let (...) denote x J VfVf exp{—Sef /}{■■■)■ 
We replace by its "dressed" version, (Fq + (C))^^, 
which includes the average effect of the conduction elec- 
trons ((C)) on the f-states themselves [Q. (C) acts like 
an f-state self-energy. Thus we write 



ln(det(F„ + C))^ 
«Trln(Fo + (C)) 



ln(det(Fo 
+ Tr((^^„ 



{C)+C)) 



(8) 



Trln(_F'o -|- (C)) is a constant which we will ignore. In 
the matrix C^^, we set Afc = and Au) — which 
correspond to SI and S2 scattering processes. The only 
nonzero components of (C) are the diagonal terms of A: 

Co = ((i?)++,o - {R)—,o)pS{Ak)6{Auj) (9) 

and the off-diagonal terms of B: 



= {R)+-^^l36{Ak)S{ALu) 



(10) 



since {R)^ — ,Tr ~ {R) — h,7r, by symmetry. 

The last term in eq. (||), together with Soc, gives our 
approximation to 5*6/ / : 



-T 



dk 



D 



lOJ 



27r 
h — 2t cos(fc)(T 

JC 



c{kuj)Dc{kuj) 



(11) 



JC„ 
4A 



4A 



■a 



— h — 2t cos(fc)(T^ 



where A = y'{h/2 - CoY+Cl and h = h/2 - J{h - 
2Co)/8A. Analytic continuation of D^^ to real time 
yields the effective conduction electron propagator whose 
poles give the excitation spectrum which consists of four 
branches. The upper two branches 



{\h\ ± |7|)2 + (C^J/2A) , where 7 = 2tcos{k), and the 
lower two branches are — e±(Co,C^). The gap between 
the upper and lower branches is A^p = Ct^J/X. Thus 
Aqp would signal a MIT. 

To determine the gap, we need C-^ and A; A is a func- 
tion of Co and C^. In eqs. (p|) and (O) Co and Ctt 



involve averages over Seff, which is itself a function of 
these parameters. As a result, we get self-consistent cou- 
pled equations for and Co- 



Co 



^ dk \ h + jsgn(h) h — ^sgn{h) 



^ dk 



£+(Co, Cjr) 
1 



^-{Co:C^) j 

^ 1 
e-|-(C'o, Ctt) e_(Co,C7r) 



(12a) 
(12b) 



Inspection of eq. (12a) shows that \Co\ < J/ 4:, which 
indicates that A will not diverge due to Co diverging. (If 
A diverges, the gap goes to zero.) To find C^, note that 
as long as the two lower bands are not completely sepa- 
rated in energy, i.e., \h\ < 2t, the right hand side of eq. 
(121) diverges logarithmically as 0. This singu- 

larity guarantees that the self-consistent equations can 
always be satisfied by a finite value of C^ for \h\ < 2t. 
Thus the gap never vanishes. (For \h\ > 2t, the gap 
grows linearly with h as discussed earlier.) This conclu- 
sion does not rely on the precise shape of the band, nor 
on TT momentum transfer for the scattering processes. 
We can solve eqs. (|12a) and (12b) numerically for the 



gap. As shown in Figs. ^ and ^ there is good agreement 
between our analytic (adapted for finite size) and DMRG 
approaches which indicates that our analytic calculation 
contains the essential physics. In Fig. ^ the analytic 
results are smaller since they include only a subset of 
all possible scattering processes. Note that the analytic 
results are for periodic boundary conditions. 

Using eqs. (^2^ and (|l2^), we find that A™" for 
an infinite lattice becomes very small as J — > (see in 
the inset of Fig. ||) |0. For a linear conduction band 
(7 = 2i|l — 2k/TT\) in the limit of small J, 



(13) 



where gi , 172 are real positive constants. The factor of 
in the exponent shows that pairs of spin flips dominate, 
e.g., an f-spin, which points up in the field, flips from up 
to down and back to up. This functional form also fits 
the gap for a cosine conduction band (7 — 2tcos(fc)) as 
shown in Fig. ^, though we could test the fit only over a 
rather limited range of J since the extremely small values 
of A™" for J < 1.8 were obscured by numerical noise. 

We argue that there is also no field induced MIT for 
the half-filled asymmetric Anderson lattice, which lacks 
particle-hole symmetry. The external field pushes the 
f-band, which has spins aligned along the field, away 
from the Fermi energy and into the Kondo regime, where 
charge fiuctuations are reduced. This can be shown using 
the Schrieffer- Wolff transformation. DMRG calculations 
on small lattices (L < 8) also do not find a MIT psj . 

We believe that spin fiip scattering will also keep the 
gap open in three dimensions. The major difference in 
3D is the formation of Landau "tubes" . We conjecture 
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that the associated van Hove singularities will actually 
lead to enhanced scattering and hence to a larger A™" 
than in the ID case. 

To summarize, ID Kondo insulators do not undergo a 
magnetic field induced MIT due to strong spin flip scat- 
tering across the Fermi surface. Using a new analytic 
approach, we find that for an infinite lattice, the quasi- 
particle gap depends exponentially on for small J. 
We suggest that the absence of a field induced MIT ap- 
plies to the 3D case as well as the asymmetric Anderson 
model. The fact that a large field can significantly reduce 
but not close the insulating gap is consistent with exper- 
iments which find that Ce3Bi4Pt3 acts like a semimetal 
in high fields ^ . 
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band when h ^ hs leads to a larger gap than when h < 
hs- (The density of states becomes large as L increases.) 
[15] Mariana Guerrero, private communication. 

FIG. 1. Quasiparticle gap versus h for J — 2 with L — 32. 
The thick solid line is the DMRG result with open boundary 
conditions, while the thin line is the analytic result with peri- 
odic boundary conditions. Inset: Typical conduction electron 
band structure for J = 0, t = 1 and h = 2. 

FIG. 2. Minimum quasiparticle gap versus J for L — S, 16, 
and 32. Solid symbols are DMRG results; open symbols are 
analytic results. Solid lines are guides for the eye. Inset: Min- 
imum quasiparticle gap versus J for an infinite lattice. The 
solid line is a fit to eq. ( p^ ) with gi = 12.13 and g2 = 57.97. 
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